Non-Markovian effects in the solar neutrino problem 
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Abstract. The solar core, because of its density and temperature, is not a weakly- 
interacting or a high-temperature plasma. Collective effects have time scales comparable to 
the average time between collisions, and the microfield distribution influences the particle 
dynamics. In this conditions ion and electron diffusion is a non-Markovian process, mem- 
ory effects are present and the equilibrium statistical distribution function differs from the 
Maxwellian one. We show that, even if the deviations from the standard velocity distribu- 
tion that are compatible with our present knowledge of the solar interior are small, they are 
sufficient to sensibly modify the sub-barrier nuclear reaction rates. The consequent changes 
of the neutrino fluxes are comparable to the flux deficits that constitute the solar neutrino 
problem. 

The solar neutrino problem is one of the most interesting long standing puz- 
zle of the modern physics. The discrepancy between the combined results from 
the solar neutrino experiments (Homestake, GALLEX, SAGE, Kamiokande and Su- 
perKamiokande) and the predictions of the standard solar models (SSM) combined to 
the minimal standard electroweak model has suggested the hypothesis that neutrinos 
could have small masses and the lepton flavor is not conserved ||]. Because the 
neutrino oscillation theory has far-reaching consequences for both particle physics and 
cosmology, it is of great importance to ask whether the solar neutrino problem can 
be solved, or at least alleviated, in the framework of the conventional physics ||. 

The aim of this contribution is to show that the physical characteristics of the 
solar core (the plasma parameter T — (Ze) 2 /(akT), where a is the interparticle av- 
erage distance, the collisional frequency v c and the electric microfield intensity and 
distribution) indicate the presence of long-range many-body interactions and memory 
effects. Therefore, the appropriate equilibrium distribution function of the ions and 
electrons deviates respect to the usual Maxwell-Boltzmann (MB) distribution, which 
is recovered in the weakly-interacting limit. These small deviations are amplified by 
the barrier penetration mechanism yielding non-negligible modifications of the solar 
neutrino fluxes 

Nuclear reaction rates in the solar interior (SI), which has a plasma parameter 
r k, 0.1, can be enhanced or suppressed depending on the behavior of the high- 
energy tail of the equilibrium distribution function of such weakly non-ideal plasma. 
Similar effects are expected in other environments characterized by plasma parameters 
0.1 <r <1, e.g., recombination and nuclear rates in brown dwarfs, Jupiter core and 
stellar atmospheres. 

In the SI the conditions required by the Debye-Huckel are only approximately 
verified. The plasma frequency v p i is of the same order of the collisional frequency v c : 
v p i ~ v c w 3 — 6 x 10 17 sec -1 and the screening radius is of the order of the interparticle 



distance. In addition to many-body collisional effects, electric microfields are present. 
The distribution and the fluctuations of these microfields must be carefully considered, 
since it modifies the usual Boltzmann kinetics. 

In the SI all these effects that can modify the Boltzmann kinetics are not very 
strong, contrary to the case of strong interacting Coulomb plasmas (r > 1) ||. 
Therefore, corrections to the usual MB distribution are small and affect mainly the 
more energetic particles. Physical quantities that depend on integrals over the entire 
distribution are almost unaffected. However, nuclear reactions are strongly affected, 
when they receive the main contribution from the high-energy tail of the distribution, 
because of the Coulomb barrier. 

In general, small deviations from the MB equilibrium distribution can always be 
asymptotically described by M 
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where 5 is the Clayton parameter Q. In particular, it has been shown 0, Q that, if 
one expresses the small corrections to the usual kinetics by a suitable power expan- 
sion of the diffusion and drift coefficients in the Fokker-Plank equation, the leading 
modification of the distribution is indeed given by Eq. (Q) . 

In a more general context, when deviations are not necessary small, it is necessary 
to generalize the Boltzmann-Gibbs statistics and to include the possibility of non- 
extensive thermostatistics. Tsallis has recently developed a statistical framework that 
naturally includes long-range interactions and non-Markovian memory effects §. The 
particle velocity distribution predicted by this theory is 



/(«) = 



I + (<?-!) 



2kT 



1/(1-9) 



e 



! + (<?-!) 



2kT 



(2) 



where G is the Heaviside step-function, and the parameter q can be connected by 
expansion to 6: q = 1 — 26. In the limit q — > 1 the Boltzmann-Gibbs statistics is 
recovered. 

A specific physical model for the microscopic processes in SI that produces the 
above equilibrium distribution is outlined in the following. 

Each particle is affected by the total electric field distribution due to the other 
charges in the plasma. The random electric microfields are often expressed in term of 
the dimensionless parameter F as (£ 2 ) = (Fe/a 2 f , where 47ra 3 /3 = l/ra; the distri- 
bution of F in plasmas depends on the value of T M [h], p] |TJ, [l3| . These microfields 
have in general long time correlations, and can generate anomalous diffusion. 

The total microfield can be decomposed in three main components, (i) A slow- 
varying (relative to the collision time) component due to the collective plasma os- 
cillations, which the particle sees as an almost constant external mean field £ over 
several collisions, (ii) A fast random component due to particles within a few Debye 
radii, whose effect can be described by an elastic diffusive cross section o\ ~ v^ 1 . 
When only this cross section is present, the distribution remains Maxwellian even 
in presence of the slow mean field £ . (iii) A short-range two-body strong Coulomb 
effective interaction, that can be described by the ion sphere model H. The strict 



enforcement of this model as implemented by Ichimaru Q yields the elastic cross sec- 
tion ao = 2ira 2 a 2 , where a is the interparticle distance, a an adimensional parameter 
whose order of magnitude can be inferred from the parameter F that characterizes 
the microfields, F » a~ 2 . Since F 2 ~ 3/r w 40, for T = 0.07, one can estimate 
0.4 < a < 1. In the present model, it is this component of the electric held that turns 
out to be mainly responsible of the correction factor exp [— 6(E/kT) 2 ]. 

In this framework, the stationary solution of the kinetic equation valid for small 
deviations from the MB distribution can be shown to be: 
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K = 2m a mb/ (rria + mf,) 2 is the elastic energy-transfer coefficient between two particles 
a and b, and S is proportional to the square of the ratio of the energy densities of the 
electric field and of the thermal motion. 

In the small correction limit, relevant to the SI, 0=1 and the Clayton parameter 

is: 

|l| ~5Rr 12alr2<<1 - <6) 

While we presented only a specific model, it is already sufficient to introduce the kind 
of mechanisms that modify the MB distribution. It is especially important to study 
the time and spatial correlation of the microfields, and how the physical phenomena of 
interest, e.g., subbarrier nuclear reactions, generate the appropriate scale separation 
of the random field contributions. 

The presence of even a tiny deviation from MB in the solar core produces large 
changes of the subbarrier nuclear reaction rates and, consequently, of the predicted 
neutrino fluxes. Following the general homology relationships for the variations of 
physical inputs, see for instance Ref. we can estimate the effect of the non-Maxwel- 
lian distribution on the fluxes || : 
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for the fluxes j = 7 Be, 8 B, 13 N and 15 0, while we use the solar luminosity constraint to 
determine the pp flux, R pp = l + 0.087x (l-i? Be ) + 0.010x (1 - i? at) + 0.009 x (1-Ro), 
and keep fixed the ratio £ = $ pe p/&pp = 2.36 x 10~ 3 . The power index (3j, that 
compares in Eq. (^), depends from the nuclear reaction considered and its value has 
been taken from Ref. In principle, there could be a different parameter Sj for each 
reaction and it should be possible to calculated them from the specific interactions 
in the solar plasma core. For the purpose of estimating the possible effects of this 
mechanism on the solar neutrino fluxes, we used two simple models where Sj were 
used as free parameters. 



In the first model, we have used the same S for all the reactions; the best fit 
to the experimental data give 6 = 0.005 with a corresponding \ 2 — 35. I n the 
second model, we have fitted two different j's, one (6(17)) for the the p + 7 Be reaction 
and the other (£(34)) for the 3 He + 4 He reaction. The best result gives % 2 = 20 
with <5(i7) = —0.018 and (5(34) = 0.030. The negative value of $(17) means that the 
corresponding distribution has an enhanced tail and that the p + 7 Be reaction rate 
increases. 

In spite of the fact that the values of \ 2 are much smaller than the ones in the 
SSM (Xssm > 74), they are still large: this mechanism cannot solve the solar neutrino 
problem. The best values of S are of the order of the ones found using Eq. (§) with 
a w 0.4 and T w 0.07, and they have been found not incompatible with the constraints 
from helioseismic observations . Even if 6 is small, it has non-trivial consequences 
on the neutrino fluxes: the boron (beryllium) flux can change of as much as 50% 
(30%). 

In conclusion we have shown that it is plausible that the equilibrium velocity 
distributions in SI are not Maxwellian and they could follow from Tsallis statistics. 
The parameter S that characterizes the generalized distribution is related to the mi- 
croscopic collective processes in the plasma. In particular, it is of great interest to 
consider the influence of the random electric microfields on the diffusion processes, 
and the possibility of characterizing their effects by average quantities (mean fields 
and cross sections) over different scales. In the future, one needs detailed studies 
of the time (memory effects) and spatial correlations of these microfields. Although 
these effects cannot by themselves solve the solar neutrino problem, they can produce 
large changes of the predicted fluxes and should be taken into account when looking 
for other solutions to the solar neutrino problem. 
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